The Harnack inequality for a class of 
degenerate elliptic operators 

Frangois Hamel a and Andrej Zlatos b 

a Aix-Marseille Universite & Institut Universitaire de France 
LATP, Faculte des Sciences et Techniques, F-13397 Marseille Cedex 20, France 
k Department of Mathematics, University of Wisconsin, Madison, WI 53706, USA 

Abstract 

We prove a Harnack inequality for distributional solutions to a type of degenerate 
elliptic PDEs in N dimensions. The differential operators in question are related to the 
Kolmogorov operator, made up of the Laplacian in the last N — 1 variables, a first-order 
term corresponding to a shear flow in the direction of the first variable, and a bounded 
measurable potential term. The first-order coefficient is a smooth function of the last 
N — 1 variables and its derivatives up to certain order do not vanish simultaneously at 
any point, making the operators in question hypoelliptic. 

1 Introduction 

We prove a Harnack inequality for distributional solutions to the degenerate elliptic PDE 

A y u + (3(y)u x + j(x,y)u = (1.1) 

in cylindrical domains in ~R N with axes in the direction of the first variable x. Here 7 is 
bounded measurable and (3 is a smooth function such that the operator 

AT-l N-l 

L = Y,K + X,:= £)ftj 2 + /W* = A v + (1-2) 

n=l n=l 

satisfies Hormander's hypoellipticity condition. That is, vector fields {X n }^=Q and their 
commutators up to certain order span the whole tangent space ~R N at each (x, y). Moreover, 
j3 changes sign so that L is not parabolic, since then the "elliptic" Harnack inequality fll.4}) 
below would not hold in general. These conditions on (3 are equivalent to hypothesis ( 11.31) 
below and our result is then as follows: 
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Theorem 1.1 Let D C M. N 1 be open connected and u : (a, 6) x D — > [0, oo) a bounded 
distributional solution of (II. ip 7 bounded measurable and (3 satisfying for some reN, 

/3 G C°°(D), inf /3 < < sup /?, and V \D c /3(y)\ > /or a// y e D% (1.3) 

0<[C|<r 

T/ien /or eac/i [a', b'] C (a, 6) and bounded open D' with D' C £), i/iere is C > 0, depending 
only on D, D' , f3 and an upper bound on (a' — a) -1 , (6 — b')" 1 , b' — a' , and H7H00; such that 

sup u < C inf u. (1.4) 

(a>,b>)xD> (a',b')xD' 

Remark. We note that A y could be replaced by any x- independent, uniformly elliptic in 
y operator on D, but for the sake of simplicity we state the theorem with A y instead. 

This result is motivated by its application in our work [6] on large amplitude A — > 00 
asymptotics of traveling fronts in the ^-direction, and their speeds, for the reaction-advection- 
diffusion equation 

v t + A a(y) v x = A x>y v + f(v) (1.5) 

on IR Ar+1 , with the first order term representing a shear flow in the ^-direction and / a non- 
negative reaction function vanishing at and 1. The front speeds in question are proved to 
satisfy lim^oo c*(Aa, f)/A = «;(«,/) for some constant «;(«,/) > 0, so after substituting 
the front ansatz v(t, x, y) = u(x — c*(Aa, f)t, y) into (II. 5p and scaling by A in the x variable, 
one formally recovers (II. ip in the limit A — > 00, with f3(y) := n(a, /) — a(y) and j(x, y) : = 
-f(u(x,y))/u(x,y). 

The study of hypoelliptic operators of the form 

M 

L = Y,X 2 n + X 

n=l 

(where X n are first order differential operators with smooth coefficients), possibly with an 
additional potential term, has been systematically pursued since Hormander's fundamental 
paper [7]. Although various regularity and maximum principle results have been obtained 
soon thereafter (see, e.g., [21 El HJ HU Q32 [2D]), Harnack inequalities and related heat kernel 
estimates for such operators have initially been proved only in the case when the tangent 
space at each point is spanned by the fields {X n }ff =1 and their commutators, sometimes with 
Xq either zero or a linear combination of {X n }^ =1 |9j HUJ H2J [13] . 

More recently, Harnack inequalities have been obtained without this assumption for 
certain special classes of operators, not including (II. ip with general (3, 7. Specifically, 
some operators with constant and linear coefficients, such as the Kolmogorov operator 
L = dy y + yd x — d t , were considered in [3 [16], and cases of more general coefficients satisfying 
somewhat rigid structural assumptions (see hypothesis [H.l] in |T7]) were studied in [TTl [T7] 
and with a potential term in [TB]. The domains involved in the obtained inequalities have to 

iFor C = (Ci, • • ■ , Gv-i) G N*- 1 , we let |C| = Ci + ■ ■ • + Ov-i and D</3(y) = (y). 
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depend on the metrics associated to the operators rather than the Euclidian metric, as shows 
a counter-example to a Harnack inequality in [5]. This is related to the need for the sign- 
changing assumption on /3 here. We also note that the operators considered in these papers 
involve the term dt and appropriate "parabolic-type" Harnack inequalities are obtained, but 
corresponding "elliptic" inequalities follow from these. 

It was a mild surprise to us that we were not able to find in the literature a sufficiently 
general result which would include our case ( II. ip . It appears that Harnack inequalities 
and heat kernel estimates become much more involved when the field X Q is required for 
Hormander's condition to be satisfied. One hint in this direction is the fact that the sign- 
changing hypothesis on j3 is necessary for (jl.4p to hold, so hypoellipticity of L is in itself not 
a sufficient condition. 

We therefore believe that our method of proof of Theorem II. II in the next section is itself 
also a valuable contribution to the problem of quantitative estimates for hypoelliptic opera- 
tors. The proof is based on the Feynman-Kac formula for the stochastic process associated 
with the operator L, and uses the independence of L, and thus also of the stochastic process, 
on x. It is not immediately obvious whether this requirement can be lifted and replaced, for 
instance, by some assumption on the relation of the stochastic processes associated to L and 
starting from two different points which can be connected by a path with tangent vector X 
at each point. We leave this as an open problem. 

We thank Luis Caffarelli, Nicola Garofalo, Nicolai Nadirashvili, Brian Street, and Daniel 
Stroock for useful discussions and pointers to references. FH is indebted to the Alexander 
von Humboldt Foundation for its support. His work was also supported by the French Agence 
Nationale de la Recherche through the project PREFERED. AZ was supported in part by 
NSF grants DMS-1113017 and DMS-1056327, and by an Alfred P. Sloan Research Fellowship. 
Part of this work was carried out during visits by FH to the Departments of Mathematics 
of the Universities of Chicago and Wisconsin and by AZ to the Faculte des Sciences et 
Techniques, Aix-Marseille Universite, the hospitality of which is gratefuly acknowledged. 



2 Proof of Theorem 1.1 



Without loss we can assume inf d> (3 < < sup D / (3 and D' connected, after possibly enlarging 
D'. We will also assume a = -5, a' = 0, b' = 1, b = 6, D = B 3 (0), D' = B^O), and H7IU < 1, 
with C then only depending on (3, because the general case is analogous. We also note that 
[TjJJ Theorem 18(c)] and boundedness of u show that u is actually continuous. 
We first claim that for each d > there is Cd,p > 1 such that 

sup u < Cd inf u, (2.6) 

[0,l]xA d ' [0,l]xBi(0) 

with A d := Aj U A~ d and A^ := {y G -Bi(O) | ± (3{y) > d}. Clearly it suffices to show this 
for all small enough d such that 7^ 0, which we shall assume. 

To this end, note that parabolic regularity theory with x as the time variable, applied on 
[-1,5] x {y e B 2 (0) I - f3(y) > d/2}, yields 

sup u <C' d o inf u, (2.7) 

[0,1] xA~ ' [2,5] xyl" 
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where C' d g > depends only on d and 0. Similarly, we obtain 

sup u <C' d a inf u, (2.8) 

[3,4] xA+ ' [-1,2] xA+ 

Next, consider the stochastic process (X x,y ,Y t x,y ) starting at (x,y) G E x ^(0) and 
satisfying the stochastic differential equation 

(dXr.dYr) = (/3(Y t X ' V )dt,V2dB t ), (X*>^Y x ' y ) = (*,v). 

Here i is a new time variable and B t is a normalized Brownian motion on M^ -1 with B = 
(defined on a probability space (f2,£>,P)). We then have 

Y^) = (X t ^ + x, >/2Bt + 2/)- (2.9) 

for any (x, y) G R x -82(0), in particular, Y*' y is independent from x. For any y G -82(0) we 
also define the stopping time 



y :=mi{t> 0\Y t x ' y £B 2 (0)}. 



T = T, 

If t A r := min{t,r}, then by the Feynman-Kac formula, || 'y II 00 < 1> an d the parabolic 
comparison principle, we have for each t > and (x, y) G K. x -82(0), 

e-*E(«TO, 1™)) < u(z, y) < e*E(«(J^, 5^ at))- (2-10) 

(The Feynman-Kac formula is usually stated for C 2 functions so we provide a proof of (I2.10p 
in Lemma 12.11 below.) Here 

E(u(X^Y% y )) = [ u(X^(u),Y^(u))dF(u) = [ u(x',y')d^ y (x',y'), (2.11) 

Jn JixB 2 (o) 



with the probability measure fi x,y onMx 8 2 (0) such that fi x ' y (A) = F((X^, Y x ^ y ) G A) for 
Borel sets KMx B 2 (0). Notice that $ y is supported on [x - \\/3\\ 00 t,x+ ||/3||oo^] x B 2 (0) 
and ^ y (R x dB 2 (0)) = P(r y < t). 

By (12. 9p . translation in x equally translates the (i x,y , and the (x- independent) measure 
on -82(0) given by u\{A) = fi x ' y (M. x A) is just the law of y/2B tATy + y, the Brownian motion 
on 82(0) starting at y, with stopping time r y , and with time scaled by a factor of two. In 
particular for each t > there is /i t > such that for any yi,y 2 G -Bi(O) and any Borel sets 
Ai C 5i(0) and A 2 C £ 2 (0), 

Mf(A 2 ) <i/f(A 2 ) <hfvf{A 2 ) and < < ^Vil- (2-12) 

From this it follows for t := H/SH^ 1 that 

inf u>C% inf u (2.13) 

[0,l]xfli(0) ' P [-l,2]xA+ 

with C% := e-*h t m.m{\A d \, \A d |} and, similarly, we obtain 

inf u>C'! B inf w. (2.14) 

[3,4]xfli(0) ' P [2,5]xA~ 



Using flUD, (EUD, d2HD, and ( 12~T3]) (in that order) yields 

sup u < Cd r inf w, 

[0,1] x A' ' [0-l]xBi(0) 

with > depending only on d and (3. An analogous argument gives 

sup u < CdB inf u, 

[0,l]xA+ ' M*^ ) 

and (12. 6p follows. 

Next we let w(x, y) := J* u(x + s, y)ds for some z G (0,1/3] 

sup v < C z g inf u (2.15) 

[0,l]xBi(0) ' [0.1]xBi(0) 

holds for some C z $ > 1. Indeed, it follows from (|279|) . (I2.10p . (12.1 ip that for each (x,y) G 
R x B 2 (0), 

e ~ l i u(x', y')dtf y > z (x', y') < v(x, y)<e l \ u(x', y')dtf> y ' z (x' , y') 



MxB 2 (0) JIxB 2 (0) 



where ^ ,y,z (x',y') = fi^ ,y (x',y') * (x[-z,z](x')dx'So(y')) ■ The above claims about /i^ and the 
definition of v\ imply that 



M 



x,y;z / i i\ ^ x:z / i\ w y i i\ ^ \ ^ a 

/V y ' (x , y ) < K t [x ) x v y t {y ) < ^ n t 



x+2rnz,y:z/ / /\ 



m=-M 



where k^' z is the measure on 1R with n^' z (B) = \BC\[x— z— ||/3||oo£, x+z + \\/3\\oot]\ for any Borel 
set B C E, and M is such that (2M+l)z > 2(z+||/?|| 0O t), for instance, M := [1/2+11/311^/^] . 
This and the first claim in (I2.12p means that 



M 

a:, 2/1) < e 2t ht 2 ^ v(x + 2mz,y 2 ) (2.16) 

m=-Al 



for any i6l, yi, y 2 G £>i(0) and t > 0. 

Now we take any x G [0, 1], y\ G i?i(0), and y 2 G for some fixed d > such that 
A± ^ 0. Pick t := (2||/3|| 0O )~ 1 ^ and M = 1 to obtain using (12~T6|) . 



2tL-2 



v(x,yi) < e M h 



/3z t-1 
u(x + s,y 2 )ds < e 2t h^ 2 / u(x\y 2 )dx'. 
-3z J-i 



Since f]2. 61) and its shifts in x give for c = — 1, 0, 1, 



sup m < Cd,/3 inf m < sup M, 

[c-l,e]xA d ' {c}xA d ' [ CiC+ i] x A d 

sup u < Cd.p inf u < Cd,/3 sup u 

[c-l,c]xA d ' {c-l}xA d ' [ c _2,c-l]xA d 



we obtain (12. 6p with [0, 1] and C^p replaced by [—1, 2] and C\g. This proves (12 .151) . Similarly, 
(EUD with [-1, 0] and [1, 2] in place of [0, 1], together with tajl . yield 



sup 

[-l,2]x_Bi(0) 



inf 

[0,l]xBi(0) 



U. 



In a similar way one can also obtain 



sup 



v < a 



inf 



u. 



(2.17) 



[-l,2]xB 2 (0) 



[0,1] x £ 2 (0) 



for some C zA(3 > (recall that B 2 (0) CC D = B 3 (0)). 

We will now need to use (jl.3p to finish the proof. This assumption makes the differential 
operator on the left-hand side of (II. ip hypoelliptic in the sense of Hormander. It follows 
that for t > 0, the measure fi x ' y is absolutely continuous when restricted to f x -82(0) and 
also to 1 x dB 2 (0) (as an (N — l)-dimensional measure in the latter case), with densities 
p t (x, y, •)» Qt(z, V, ■, ■) > such that 

Pt(z, y, x', y') = p t (0, y, x - x, y'), 

Qt(x, y, x', y') = q t (0, y, x' - x, y'), 

and pt,qt are bounded functions when restricted to y G -Bi(O) (with y' G -82(0) for p t and 
y' G 9^2(0) for q t ). For p t this follows from the same claim for the corresponding measure 
Jl^' y on M N given by (12.111) with t in place of t A r and /3 smoothly extended to a periodic 
function on M^" 1 (whose density is smooth in all arguments, [H Theorem 3]). This is because 
Jl x t ' y (A) > fi x t ' y (A) for any Borel set A C R x 5 2 (0). 

For this would follow from the same claim for the corresponding escape measure jl* ,y 
on R x dB 2 (0) given by (12. lip with r = r y in place of i A r. We know of such a result for 
bounded domains only [H Corollary 2.11] but since fi x,y is supported on a bounded cylinder, 
it applies in our case as well. Specifically, take any a_ < — ||/3||t»£ an d a+ > \\ft\\oot- There is 
a convex open domain G with a smooth boundary whose intersection with [a_,a + ] x R^ -1 
is [a_,a + ] x B 2 (0), and the intersection with [(— 00, a_) U (a + , 00)] x R^ 1 are two smooth 
"slanted" conical caps G± C R x S 2 (0) over the (N — l)-dimensional balls {a±} x -82(0) 
with the two (rounded) tips at points with y' coordinates y' ± such that ±/3(y'±) > and 
sufficiently long so that for any (x',y') G dG± D dG, the unit outer normal vector n(x',y') 
to dG± at (x', y') satisfies 



Then G satisfies the hypotheses of [H Corollary 2.11] (it satisfies the escape condition and all 
points of dG are r'-regular). It follows that the escape measure ffi y has a density q T (x, y, •, •) 
which is a continuous function of (x, y, x', y') G Gxd*G, where d*G is the set of "good" points 
of dG, that is, all (x', y') G dG except of the two cone tips, where n(x', y') = (±1, 0, ■ ■ • , 0). 
Thus q T is bounded on S := {0} x B\ (0) x (a_, a + ) x 8B 2 (0). Since {X®' y } s < tAr almost surely 
stays in (a_,a + ), we obtain g 4 < q T on 5 and g ; = on [{0} x 8x(0) x R x dB 2 (0)] \ S. 
Finally, qt(x, y, x', y') = q t (0, y, x' — x, y') shows that q t is bounded on R x -Bi(O) x R x dB 2 (0). 



n(x',y')- (!,(),■■ ■ 



0)| < KII/3H- 1 + 1) whenever ± f3{y') < 0. 
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Let d > be such that ^ 0, let z := 1/3, t := and 

Ct := max{ sup p t , sup q t } < oo. 

RxBi(0)xRxB2(0) Rx_Bi(0)xRx9_B 2 (0) 

Then pt(x, y, x', y'), qt(x, y, x', y') < C t X[x-i,x+i](x') because the measure (j% ,v is supported on 
[x - 1, x + 1] x B 2 (0), so we obtain from ([2~TU]) and (EOT]) 

sup u < C t e l I u(x' ,y')dx'dy' + C t e t / u(x' ,y')dx'dy' 

[0,l]x-Bi(0) i[-l,2]x_B 2 (0) V[-l,2]x9B 2 (0) 

< 10C t e* sup u 

[-l,2]xB 2 (0) 

< lOCtC^e* inf it 

[0,l]xBi(0) 

by using [-1, 2] = [-1, -1/3] U [-1/3, 1/3] U [1/3, 1] U [1, 5/3] U [4/3, 2] and fl2TT7j) . This is 
(11.4)1 . so the theorem will be proved once we establish ( 12.101) . 

Lemma 2.1 // u, /3,j, Xt' y ,Y* ,V ,r y are as in the proof of Theorem li.il particular, 
INIoc < U f/ien (I2TT0D holds for (x,y) GRx B 2 (0). 

Proof. Let Z x ' y = t so that = di and K := A y + f3(y)d x + <9 2 is the generator of the 

process {X x t A ', Y t x ' v , Z x ' y ). If we let v(x, y, z) := e z u(x, y), then Kv > on 1 x 5 3 (0) x R in 
the sense of distributions, that is, 




vK*(f> dxdydz > 



for any G C °°(1R x B 3 (0) x E), with K* := A y - f3(y)d x - d z the adjoint of K. 

For any e > let 5 £ G (0, l/2y/N - 1) be such that |/%) - /%')| < £ 2 whenever 
y, y' G -85/2(0) and |y — y'\ < y/N — 1 S £ . Let 77 : R — » [0, 1] be a smooth non-negative 

function supported in [—1,1], with J_ 1 r](x')dx' = 1 and ||t/IIoo < 2. For e > define the 
mollifier 

N-l 

rj £ (x, y, z) := e~ 2 5 l - N r] (-} rj (-) JJ rj 

n=l 

and let v e := v * if and <p £ ' ,x,v ' z (x' , y', z') := r] £ (x — x',y — y', z — z'). For e G (0, 1) the smooth 
function v 6 then satisfies on K x ^(0) x K 

(i^ E ) (x, y,z)= [ vK*<j> £, ' x ' v > z dx'dy'dz' 

JRxB 3 (0)xR 

+ f v(x', y', z')W) - P(y)}<f> £ x ?> y ' z (x', y', z') dx'dy'dz'. 

«xB 3 (0)xM 

The first integral is non- negative. The integrand in the second vanishes when [a/ — x\ > e or 

Wn ~ Vn\ > &e f° r some n or \z' — z\ > e, and |0^.f ' y ' z (x', y', z')\ < 2e~ 3 5l~ N , so we have 

(Aru e )(x,y,z) > -2 Ar+2 ee z+e ||n|| 0O . 
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We next apply Dynkin's formula [15j Theorem 7.4.1] to the smooth function v e , the 
process (X^' v , Y t x,v , Z^' y ), and stopping time t A r (with r = r y ), to obtain 



E 1™, «)] = 0) + E 



iAr 



Since V s — y v uniformly on [x — \\/3\\oot,x + ||/3||oo^] x -82(0) x [0,t] as e — y (by continuity 



of v) and Z^ T < t, it follows that 



e'E [n(X^ r , ^)] > E [v(X^ Ytf, Z™)] > u(x, y). 

This is the second inequality in ( 12 .101) . The first inequality is obtained in the same way, this 
time with v(x, y, z) := e~ z u(x, y), so that Kv < on R x -83(0) x M and 



(Kv e )(x,y,z) < 2 N+1 ee 



N+2 r .-z+e\ 



U 



□ 
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